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(IMPUTATIONS  OF  JET  IMPINGEMENT  ON  A FLAT  SURFACE  t 
A.  Rubel* 

Grumman  Aerospace  Corporation 
Bethpage,  New  York  11714 


Abstract 

An  incompressible,  inviscid  rotational  flow  model 
is  used  to  represent  the  normal  impingement  of  axisym- 
metric  Jets  and  the  oblique  impingement  of  two- 
dimensional  jets  upon  a flat  surface.  The  Poisson- 
type  equation  that  results  from  this  formulation  is  cast 
into  suitable  finite  difference  form  and  solved  by  relaxa- 
tion techniques.  Results,  in  terms  of  the  ground  plane 
pressure  distribution,  centerline  velocity  decay  and 
oblique  impingement  flow  field  structure,  compare  well 
with  observations. 

Introduction 

During  lift-off  of  VTOL  aircraft,  turbulent  lifting 
jets  mix  with  the  surroundings  and  impinge  upon  the 
ground,  producing  wall  jets  which  interact  to  form  up- 
wash  regions  of  high  static  pressure.  These  phenomena 
and  their  induced  secondary  flows  produce  lifting  forces 
on  the  aircraft  as  well  as  suckdown  and  attendent  lift 
losses.  An  understanding  of  the  behavior  of  these  com- 
plex flow  fields  is  required  for  the  optimization  of  \’TOL 
design. 

Efforts  to  analyze  the  overall  nature  of  the  inter- 
actions via  the  numerical  solution  of  the  Navier-Stokes 
equations  and  a single  equation  turbulence  model  have 
achieved  some  success^'^  but  there  are  difficulties 
associated  with  the  specification  of  a unified  length 
scale  over  the  entire  computational  domain.  Incorpora- 
tion of  higher-order  turbulence  models  and  computa- 
tional resolution  of  the  three-dimensional  aspects  of  the 
flow  field  require  considerably  more  time-  consuming 
calculations  and  have  not  yet  been  attempted. 

An  alternative,  component  approach,®"®  proceeds 
by  dividing  the  flow  field  into  separate  interacting 
regions  which  are  coupled  to  provide  an  overall  solu- 
tion. The  advantage  of  this  approach  is  that  each 
region  may  be  represented  by  models  closely  related 
to  the  physics  of  the  problem  and  relatively  accurate 
engineering  solutions  obtained  with  comparatively 
simple  methods.  Within  the  context  of  Ais  approach, 
scaling  analyses®’®  indicate  that  the  details  of  the 
impingement  region  are  dominated  by  a balance  of  the 
pressure  forces  due  to  flow  deflection  and  the  inertial 
forces  of  convection. 

In  this  study  a variety  of  incompressible  impinge- 
ment problems  are  considered  under  the  assumption 
that  the  flow  within  the  impingement  zone  can  be  treated 
as  an  inviscid,  rotational  flow.  Axisymmetric  jets 
impinging  normally  and  two-dimensional  jets  impinging 
obliquely  form  a set  of  problems  that  can  be  formulated 
in  two  dimensions  with  vorticity /radius  or  vortlcity, 
respectively,  conserved  along  streamlines.  Both  de- 
veloping and  fully  developed  jet  velocity  profiles  are 
considered  as  entrance  conditions  to  the  impingement 
region  and  finite  difference  techniques  are  used  to 
solve  the  governing  equations.  The  results  are  compared 


with  observations  in  order  to  determine  the  range  of 
validity  of  the  inviscid  rotational  model  with  respect  to 
the  pr^iction  of  the  ground  plane  pressure  distribution 
and  the  structure  of  the  oblique  impingement  flow  field. 


Formulation 

The  normal  impingement  of  axisymmetric  jets  and  the 
normal ''oblique  impingement  of  two-dimensional  jets  are 
represented  here  by  the  equations  governing  inviscid 
rotational  flow.  The  jets  arc  considered  to  be  far  from 
the  Impingement  plane  so  that  influx  conditions  are  not 
affected  by  the  jet  deflection.  The  formulation  is  given 
in  terms  of  dimensionless  quantities;  velocities  are 
scaled  by  the  maximum  jet  velocib  , distance  by  the  jet 
half  width  (i.e. , distance  measured  along  the  perpen- 
dicular from  the  jet  centerline  to  the  point  at  which 
the  velocity  is  half  the  maximum)  and  pressures  by  the 
total  pressure  at  the  jet  centerline.  Ambient  and  jet 
static  pressure  are  taken  to  be  zero  without  loss  of 
generality. 

A coordinate  system  is  chosen  with  z directed  per- 
pendicular to  the  ground  plane,  positive  outward,  and 
X directed  along  the  ground  plane  such  that  the  origin 
is  at  the  intersection  of  the  ground  plane  and  the  jet 
centerline  (Fig.  1).  Eliminating  pressure  from  fte 
X and  z momentum  equations  ylel^  the  vorticity 
equation 
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Fig.  1 Schematic  of  the  Jet  Impingement  Problem 
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d(t«>/x°)  = 3^^°)  + w = 0 (1) 

dt  dx  Bz 

where  u and  w are  the  velocities  io  the  x and  z direc- 
tions, respectively,  n is  an  index  of  the  flow  character 
(axisymmetric  n = 1,  two-dimensional  n = 0)  and  cj  is 
the  vorticity  defined  by 

It  - If  • <2) 

A dimensionless  stream  function 

8 1/'  a a 

assures  satisfaction  of  mass  conservation  and  when 
combinc<l  with  the  definition  of  vorticity  produces  the 
Poisson-type  equation  that  provides  the  basis  for  the 
ensuing  analysis 
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The  unifying  feature  of  the  flows  considered  here  is 
that  the  vorticity  function,  tj  'x" , is  conserved  along 
stream  surfaces  and,  in  the  absence  of  separation 
buljbles,  is  fully  specified  in  terms  of  the  incoming 
jet  stream  function.  The  pressure  is  determined 
from  an  auxiliary-  Bernoulli  equation 
■>  ■> 

p • u“  * w - f ( 1 ( 


which  follows  from  the  vector  momentum  equation. 
This  relationship  implies  that  the  total  pressure  is 
conserved  on  stream  surfaces  and  completes  the 
formulation  of  the  problem. 


Although  a unified  treatment  of  the  class  of  flows  to 
be  investigated  has  been  presented  to  this  point,  it  is 
adv;mtageous  to  further  categorize  the  flows  to  more 
clearly  indicate  the  manner  in  which  boundary  condi- 
tions are  prescribed  and  special  cases  are  developed. 
For  normal  impingement  there  is  flow  symmetry  about 
the  jet  i-cnterline  (x  = 0)  and  only  the  domain  x ■ 0 
need  be  i-onsidered.  The  boundary  conditions  are 


4 (0,  z)  = 0 

symmetry 

(6  a) 

(oc,  Z)  = 0 

parallel  outflow 

(Gb) 

^(x,  0)  = 0 

ground  plane 

(6c) 

4{x,oc)  = 4{\) 

jet  inflow 

(6d) 

The  boundary  condition  at  x—  «o  can  be  replaced 
by  a stream  function  distribution  ♦(z)  (e.g. , Para- 
meswaren^O)  and  it  is  Instructive  to  demonstrate  the 
relationship  between  this  distribution  and  condition  (6d) 
as  well  as  the  implications  of  condition  (6b). 

Consider  the  points  (X,  z)  for  X -•  » and  (x,  Z)  for 
Z ->  « in  order  to  determine  a relationship  between 
X and  z such  that  these  points  lie  on  the  same  stream- 
line. The  boundary  conditions  (6b,  d)  imply  that  when 
X—  •«  (Z-*  «)  the  second  (first)  term  of  equation  (4) 
becomes  negligible  and  x and  z are  related  para- 
metrically through  w/x“(^>).  For  the  axisymmetric 
case,  then 


x^Lax^Vaz;  ax  3,2j  ^23^2  ^3  ax 

and,  equating  coefficients  of  the  stream  function 
derivatives, 

X .®JL  = A =-Si, 

* ,2  az  V dz)  \azj  2 


3x  _ X 
az  X 

which  has  as  its  solution  x^  = 2Xz.  Thus,  for  the 
axisymmetric  case, 

fiX,  z)  = Il>(  J~25G..  Z)  (7a) 

and  in  a similar  fashion  it  can  be  shown  that  for 
two-dimensional  normal  impingement 

V-IX,  z)  = (/>(z.  Z).  (7b) 

The  boundary  condition  (6b)  can  be  recovered  by  noting 
that 

= 0 

\ax/  ..  for  the  two-dimensional  case  and 

(II)  - -zw(.rw,  z) 


for  the  axisymmetric  case.  For  the  velocity  profile 
considered  here  (Table  1)  it  can  be  seen  that 


(^\  . 5 - ^ _L_ 

\ax/x_  «,  (1+2c2Xz)^  ^*^2^ 


so  that  the  par.allel  flow  boundary  condition  is 
compatible  with  specification  of  the  far  field  velocity- 
profile. 


Several  jet  profiles  are  chosen  (Fig.  2)  in  order  to 
evaluate  the  effect  of  fully  developed  vis  a vis  develop- 
ing flow  on  the  ground  plane  pressure  distribution.  The 
functional  form  of  these  profiles^^  and  the  derived 
stream  function  distributions  and  vorticity  functions  are 
given  in  Table  1.  The  vorticity  forcing  function  is 
implicitly  related  to  the  stream  function  in  the  two- 
dimensional  developing  jet  cases,  hence,  recasting 
equation  (4)  in  a new  dependent  variable  t(^),  defined 
in  Table  1,  is  desirable.  For  these  cases  the  govern- 
ing equation  becomes 

ax2  az^  (l^erft)  LUx;/  U-/  ^ 

where  a is  a spreading  parameter  for  the  profiles  and 
the  boundary  conditions  are 

t(o,  z)  = ^ (9a) 

H (->.  z)  = 0 (9b) 

t(x,  o)  •=  (9c) 

{(X,  «)  - (l-x)/ff  (9d) 
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Fig.  2 .let  Velocity  Profiles 


Equation  (4)  remains  valid  for  the  oblique  two- 
dimensional  impingement  cases  (only  fully  developed 
profiles  are  considered)  but  the  independent  vari^les 
are  transformed  so  that  the  new  z coordinate  is 
aligned  with  the  Jet  centerline.  The  new  x coordinate 
is  measured  from  the  jet  centerline  in  a direction 
parallel  to  the  ground  plane.  Therefore, 

z = z/sln  /?  , X = x-z/tan3  (10) 

and  equation  (4)  becomes 


9^*  •>  a 

;72  - 

a Z 


2Cj  sin^/S<(’(l-Cj<'^)  = 0 
(11) 


where  p is  the  angle  between  the  jet  centerline  and  the 
ground  plane  (Fig.  1).  The  stream  function  is  still 
related  to  the  veloeitics  u and  w by  equation  (3)  and, 
for  0 - 90',  equation  (11)  reduces  to  equation  (4)  with 
X = X,  z = z.  The  bouixlary  conditions  for  the  oblique 
impingement  problem  are 


a ii> 

^ (±  » . J)  0 (12a, b) 

iKx,  «>  ) ■-  (CjjcsintJ)  (12c) 


(12d) 


where  ifg  is  the  value  of  the  stagnation  stream  function. 
It  is  easily  demonstrated  that,  for  the  oblique  impinge- 
ment case. 


(KX,z)  =(('  (xg  + z,  Z)X-«,  (13) 

where  Xg  is  the  distance,  in  the  x coordinate,  from  the 
jet  centerline  to  the  stagnation  streamline  at  the  jet 
exit.  Using  this  description  of  the  ground  plane  flow 
far  from  the  origin,  a balance  of  momentum  in  the  x 
direction  yields  the  value  for  t/ig.  Thus, 

2cos^.f'^'“Wd  !('=  /^”nd((i+ 

where  W and  II  are  the  velocities  in  the  jet  centerline 
(influx)  and  ground  plane  (efflux)  directions,  respectively. 
Utilizing  the  fully  developed  profile  information 

i.e.  ,-W  = ir=l-c^(f^,  tp  = 
produces  the  desired  result 

cos  13  =1  Cj«g  (l-ci[i('2/3)  (14) 

and  completes  the  formulation  for  the  impingement 
problems  considered  here. 

Method  of  Solution 

Finite  difference  techniques***  are  used  to  solve  the 
governing  Poisson-type  equation  for  the  stream  function 
field.  Since  the  decaying  influx  and  efflux  velocity  pro- 
files far  from  the  origin  are  infinite  in  extent,  the  calcu- 
lation of  a bounding  jet  streamline  is  not  of  concern. 
Instead,  the  extent  of  the  physical  domain  chosen  for  the 
application  of  the  boundary  conditions  must  be  large 
enough  that  the  results  in  the  ground  plane  are  con- 
verged (i.e. , a larger  domain  produces  no  change  in 
results).  The  normal  impingement  problem  is  solved 
in  the  rectangular  region  of  height,  Z,  and  width,  X 
(Fig.  3a).  Furthermore,  to  match  the  stream  function 
at  the  comer  point  (X,  Z),  requires  that  X = 2Z  for 
the  axisymmetrlc  case  (Equation  7a)  and  X = Z for  the 
two-dimensional  case  (Equation  7b).  (Numerical  ex- 
periments indicate  that  such  a domain  produces  improved 
mesh  conver(^nce  properties.)  The  oblique  impinge- 
ment problem  is  solved  in  the  numerically  rectangular 
region  of  height,  Z,  and  width,  23i  (Fig.  3b).  The 
displacement  of  the  stagnation  streamline,  Xg,  is  small 
enough  in  the  problems  considered  here  that  at  the 
corner  points  (±X,  ?)  the  stream  function  is  nearly 
matched  (Equation  13). 


Table  1 Rotational  Jet  Flow  Parameters 


Velocity 

-w(x) 

Stream  Function 
ij-(x) 

Vortlclty  Function 

Axisymmetrlc 

Fully  Developed 

2 -2 

(1  + c x-") 

C2= 

' — — -i 

x^/2(l  + Cg  x^) 

4c2  (1-2c2i(()^ 

Two-  Di  menslonal 

Fully  Developed 

(1  - tanh^  CjX) 

Cj  =■  tanh"*  (l/^/2) 

— tanh  (c,x) 

Cl  ' 1 ' 

2cf  *(l-ci  ^2) 

Two-Dimensional 

Developing 

j-  1 l+erf[(l-x)/ol| 

0 » ,1,  .25,  .5 

,(1-X)/<T 

■f/  (l+erff)d4 

t=  (1  - X)/<T 

0-^ 
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The  numerical  representation  of  equations  (4)  and  (8) 
makes  use  of  a uniform  mesh  with  spacing  defined  hy 


A similar  representation  of  equation  (11),  for  oblique 
Impingement,  makes  use  of  a uniform  mesh  with  spaclngs 


“ N*-l  N^-l 


AX  = 


2X 


, Az  = 


Ng-1 


where  Nx  xnd  are  the  number  of  points  in  the  x and  z 
directions,  respectively.  Using  central  differencing 
for  the  derivatives  of  the  dependent  variable  results  in 
the  internal  point  (1,  j)  finite  difference  equations  for  the 
axisymmetric  case  (l.e. , Equation  4,  n = 1), 


- <^1.1,)*  12*“"  '’'-i/2"  V2' 

-Sc^xf  j (1-2c2«,  j)  j 


* 402^^  I (l-'Jc^^l'j^  j)  (Az)‘, 

where  I)  ^ Az  'Ax  and  x , = x.  ,/(x.  . i Ax/2),  for 
) 2 i.J  ' 1.) 

the  two-dimensional  fully  developed  case  (l.e.,  Equation 
4,  n = 0), 


- <^'.1^,- »-  (♦l.j.l’  '/'i.j-l'  (l(i) 


and  for  the  two-dimensional  developing  flow  case  (l.e. , 
Equation  8), 


" '^i.J-1  " ^i,j-l*  * ‘’l.j[T“'^l,)'l  ‘^i.j-l* 

- (Az^o)"]  (17) 


U . 

where  A,  . = -j^(«l,l  j - «,.i  ,). 


-« 


13  = 

iJt  j mn 


1.) 


l.j  ^/^  (1  * erf  Cjj)  ' 


These  equations  are  modified  slightly  at  x = X by  the  use 
of  additional  fictitious  |x>ints  to  satisfy  the  gradient 
boundary  condition  by  reflection.  The  equations  (15-17) 
are  arranged  for  line  relaxation  by  tridiagonal  solution 
techniques.  The  asterisk  superscript  refers  to  an 
intermediate  value  of  the  dependent  variable;  coefficients 
of  these  variables  and  the  right  sides  of  the  equations  refer 
to  previous  values  of  these  quantities  (i.e. , iterative 
values).  The  (v+l)  values  are  obtained  by  the  relaxation 
formula  (18) 


where  0<  r < 1 indicates  under- relaxation  and  1 < r < 2 
indicates  over-relaxation.  ITie  nonlinear  vorticlty  forcing 
functions  have  been  arranged  in  equations  (15-17)  so  that 
they  are  partially  included  in  the  implicit  line  relaxation 
calculation. 


where  N2X  N2  are  the  number  of  points  in  the  5?  and  Z 
directions,  respectively  and  yields  the  internal  point 
finite  difference  equation 


2 [i . dV‘ - 

““»»•>“]  •m  - Cl.J 


(18) 


° <'''i,hi"\j-i>-T  “®^"''1+I,jn"'^i-I,j-1 


■'^i+i,j-i  ■‘^1-1, j-i) 
where  D = Az/Ax. 


The  ground  plane  velocity  is  a derived  quantity  of 
considerable  importance  since  it  determines  the  ground 
plane  pressure  by  Bernoulli's  equation.  From  equation 
(3),  and  the  appropriate  transformations,  this  velocity  is 
given,  to  second  order  in  mesh  height,  by 

2 

u(o,x)=[i(i(x,Az)  -^^(x,©)  (Az)^] /x"az  (19a) 

d z 

for  fully  developed  normal  impingement,  by 

u(o,  X)  = -|•(1  + erf  o')  [1  -<r|(x,Az)]  /Az  (19b) 


for  the  normal  impingement  developing  flow  and  by 


u(o,x)  = - [(/'(X,Az)  - 


(19c) 


for  two-dimensional  oblique  impingement.  The  value 
of  the  second  derivative  of  the  stream  function  can  be 
found  from  the  governing  differential  equations  evaluated 
at  the  ground  plane.  For  the  fully  developed  normal 
impingement  cases  this  is  easily  done 

0 in  2D  flow,  = Xw  in  axisymmetric  flow) . 

8z^  az^  / 

The  oblique  impingement  case,  however,  requires 
calculation  of  a row  of  cowled  fictitious  points  located 
at  z = - Az  (via  Equation  fS)  in  order  to  represent  the 
stream  function  second  derivative. 


(X,  o)  = |((i(x,Az)  - 2((i 

az“  ® 


l|>(x,  - Az)|/(Az)^  . 

The  ground  plane  pressure  is  evaluated  from 

p(x,o)  = u^  (x^,  Z)  - u2(>t,o).  (20) 

A considerable  number  of  calculations  are  necessary 
(e.  g. , Table  2)  in  order  to  determine  the  character  of 
the  numerical  solution.  Initial  estimates  for  all 
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Table  2 Numerical  parameters  for  selected  impingement  computations. 


Axisymmetric 

Two 

-Dimensional 

Normal 

Normal 

Oblique  (Fully  Developed) 

Fully 

Fully 

a 

0 (deg) 

Developed 

Developed 

.5 

.25 

.1 

80 

70 

60 

50 

40 

30 

X(2  5f)* 

8 

5 

5 

5 

5 

10 

10 

12 

12 

10 

20 

Z(Z) 

4 

5 

5 

5 

5 

5 

5 

6 

6 

8 

10 

81 

41 

41 

41 

41 

41 

41 

81 

41 

41 

41 

Nj,(Ny) 

101 

41 

41 

41 

41 

21 

21 

41 

21 

21 

21 

€ 

lO'*' 

io-« 

10-s 

10-s 

10-5 

10-5 

10-5 

10-5 

10-5 

10-5 

10-5 

U 

438 

370 

338 

428 

1938 

741 

816 

3401 

1003 

1435 

1799 

I 

t 1 ’Quantities  in  parentheses  refer  to  the  oblique  impingement  caes 

' ' X = width,  7.  - height,  N = number  of  points,  € = tolerance,  v = number  of  iterations 
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computations  reported  here  are  obtained  by  iinear 
interpolation  between  the  ground  and  jet  exit  stream 
function  conditions.  The  successive  relaxation  procedure 
(Equation  18)  is  applied  until  a convergence  given  by 


1.  j ^1.  j 


i,  j ' max 


(21) 


is  obtained.  It  is  found  that  the  convergence  parameter 
e = 10~(>  is  sufficient  so  that  when  c is  reduced  by  an 
order  of  magnitude  the  maximum  velocity  variation  at 
a ground  plane  point  is  less  than  10~^. 


The  domains  of  integration  and  mesh  widths  indicated 
in  Table  2 for  fully  developed  normal  impingement 
reproduce  the  parallel  flow  conditions  at  x = X (i.e. , 
Equations  (7a,  b)  arc  satisfied).  It  can  be  shown  that 
the  criterion  (X.lz)^'^ ' 1 must  be  satisfied  in  addition  to 
condition  (7a)  in  order  for  the  axisymmetric  ground 
plane  velocity  tft  achkwe  its  asymptotic  value,  u(X,0)— 1. 
The  case  tabulated,  for  example,  has  a slight  overshoot 
in  velocity  Ixiginning  at  x = 3,7  and  reaching  a maximum 
of  1.6%  at  X = X = 8. 


As  might  be  expected  from  a glance  at  Fig.  2,  the 
developing  flow  cases  are  the  most  difficult  to  compute. 
The  resolution  of  the  Jet  velocity  profiles  with  a uniform 
mesh  becomes  impractical  for  spreading  parameters, 

O’  s 0.1.  Furthermore,  the  nature  of  the  variable,  | , 
is  such  that  it  grows  negatively  without  bound  as  the 
width  of  the  numerical  domain  is  increased.  A cut-off, 
such  that  4 2 - 3.49,  is  used  to  restrict  this  growth  and 
only  points  satisfying  this  criterion  are  considered  within 
the  domain  governed  by  equation  (17).  This  procedure, 
along  with  appreciable  under-relaxation  (e.  g. , r = 0, 3 
for  0=0.1),  provides  converged  solutions  which  recover 
the  asymptotic  ground  piane  velocity  even  for  the 
smallest  spreading  parameter  (i.e.,  u(0,X)  = 0.98  at 
o=  0.1). 


Convergence  in  the  oblique  impingement  case 
requires  a domain  that  increases  as  the  inclination 
angle,  0,  departs  from  90*.  Numerical  experiments 
indicate  that  2 = 5/sin  0 and  X = Z (from  Equation  13) 
provides  an  adequate  domain  for  these  computations. 
The  mesh  width  convergence  at  0 = 60*  is  such  that  a 
41  X 21  grid  produces  ground  plane  velocities  within 
1%  of  those  predicted  by  an  81  x 41  grid.  A similar 


result  is  found  for  the  normal  impingement  case.  In 
order  to  further  estimate  the  convergence  properties  of 
the  oblique  impingement  cases  of  Table  2,  the  resultant 
force  on  the  ground  plane  and  the  moment  about  the 
origin  are  computed  from  the  numerical  distribution  of 
the  ground  plane  pressure  and  compared  with  their 
exact  values  assuming  recovery  of  the  initial  profiles  at 
large  X (i.e. , Equation  13  is  satisfied).  Thus, 

1 ^ 

F (calc)  = Y / jj  P 
F (exact)  = - 2 Bln0  j "wd^/i 


M (calc)  =-5-/  pxdx, 

^ -X  (23) 

M^(exact)  = Xg  F (exact). 


From  a moment  balance  about  xg,  the  location  of  the  re- 
sultant force,  it  can  be  seen  by  inspection  that  xs  is  also 
the  distance  from  the  Jet  centerline  to  the  stagnation 
streamline  measured  at  the  Jet  exit  and  along  the  direc- 
tion parallel  to  the  ground  plane  (Fig.  1). 

It  is  demonstrated  (Table  3)  that  the  computed 
force  and  moment  are  in  error  by  at  most  3%  and 
2%,  respectively,  indicating  that  the  oblique 
Impingement  cases  are  sufficiently  converged. 

The  computations  require  about  10”^  seconds/lter- 
ation/(num):«r  of  points)  and  have  been  performed  using 
an  HP  3000  Series  n computing  system  (a  factor  of 
about  20  relates  this  time  to  IBM  370-168  time).  Few 
experiments  have  been  done  to  optimize  the  relaxation 
parameter;  the  cases  of  Table  2,  except  for  the  develop- 
ing profiles,  used  r = 1.  The  cases  listed  in  Table  2 are 
those  used  to  demonstrate  convergence  properties  and 
quite  satisfactory  solutions  can  be  obtained  for  many  of 
the  cases  with  considerably  less  stringent  conditions. 
Finally,  although  the  cases  presented  here  include  im- 
pingement angles  down  to  30*,  the  limitation  is  on  pl^si- 
cal  rather  than  numerical  grounds;  converged  computa- 
tions have  been  established  for  0 as  low  as  10*. 
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Tabic  3 Comparison  of  computed  and  exact  ground 
plane  force  and  moment  as  a function  of 
impingement  angle. 


0 

F (calc) 

F (exact) 

Mo  (calc) 

Mo  (e.xact) 

80’ 

1.  .■>35 

1.490 

-. 1972 

-.2005 

70 

1.45ii 

1.422 

-. 3958 

-.4059 

liO 

1.  121 

1.310 

-.0120 

-.0219 

50^ 

1.  170 

1.  1.59 

-.8.100 

-.8503 

40’ 

0.977 

0.  972 

-1.  140 

-1. 121 

30 

0.  750 

0.750 

-1.402 

-1.430 

a)  Normal  Impirtgamant 


in  the  undisturbed  free  jet  at  the  identical  distance  from 
the  nozzle  exit.  According  to  the  present  computations 
the  location  of  this  point  is  at  z » 2. 1 which  is  in  close 
agreement  with  the  value,  z =«  2.2  observed  by  Giralt 
et  al.  Comparison  of  the  calculated  ground  plane 
pressure  distribution  with  data^>^^  (Fig.  4)  shows 
excellent  agreement  over  the  entire  range  of  the 
distribution.  This  is  a considerable  improvement  upon 
the  inviscid  solution  of  Scholtz  and  Trass*-^  which 
breaks  down  at  large  distances  from  the  centerline. 

The  agreement  between  the  current  computations  and 
those  ol  Scholtz  and  Trass  for  velocity  decay  along  the 
centerline  is  demonstrated  in  Fig.  5:  the  discrepancy 
with  respect  to  the  data*^^'  I'*  is  attributed  to  mixing 
phenomena  within  the  impingement  zone*^2_ 

Comparison  of  the  results  for  the  normal  impinge- 
ment of  fully  developed  two-dimensional  jets  proceeds 
in  a similar  manner.  Numerical  computations  yield 
the  98'i,  criterion  at  z ■=<  3.40  which,  in  conjunction 
with  accepted  spreading  ratcs^l,  provides  the  relation- 
ship  h =■  12.2,  where,  again,  h is  tlie  dimensionless 
distance  from  the  jet  exit  to  the  ground  plane.  The 
calculations,  (Fig.  0),  indicate  tliat  for  this  case,  too, 
there  is  excellent  agreement  with  observations. 

If  the  jet  exit  plane  is  in  proximity  to  the  ground 
plane  (but  outside  the  impingement  zone),  then  a core 
region  is  present  and  a developing  velocity  profile 


Fig.  3 Physical  Domain  and  Houndarv  Conditions 
for  Finite  Difference  Analysis 


Results  and  Discussion 


The  computations  described  here  are  for  idealized 
inviscid  rotational  jets  located  far  from  the  impinge- 
ment plane  with  lengths  scaled  by  the  jet  half  width  and 
pressures  by  the  maximum  jet  total  pressure.  Heal 
Jets,  however,  encounter  a free  mixing  zone  prior  to 
entry  into  the  impingement  region  and  data  are  usually 
presented  with  length  scales  given  b>'  cither  the  jet 
height  above  ground  or  the  half  width  of  the  undisturlted 
jet  at  the  ground  plane.  If  a matching  point  iK'tween 
the  computed  solution  and  point  of  free  jet  entry  to  the 
impingement  zone  is  established  then  the  free  jet  length 
scales  can  lx-  related  to  the  computational  length  scale 
liy  well  known  free  jet  spreading  rate  parameters®>  ^2, 13, 
In  presenting  the  computation/data  comparisons  it  is 
convenient  to  use  the  rotational  jet  length  scale  so  that 
dimensionless  jet  heights,  h,  arc  always  referred  to 
this  quantity.  The  ground  pressure  data  arc  normalized 
Ijy  the  peak  ground  plane  pressure,  p^,  to  eliminate 
total  pressure  loss  consi^rations. 

For  normal  oxisymmetric  impingement,  Uiralt  cl  al. 
use  the  criterion  that  the  impingement  region 
commences  when  the  centerline  vclocitj'  is  98'o  of  th.at 


Fig.  4 Ground  Plane  I’ressure  Distribution, 
Axlsymmetric  Normal  Impingement 
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Fig.  5 Centerline  Velocity  Decay, 

Axlsymmetric  Normal  Impingement 
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exists  at  the  beginning  of  the  impingement  zone.  The 
computations  of  ground  plane  pressures  under  such 
conditions  are  shown  in  Fig.  7 and  are  seen  to  fall 
between  the  potential  flow  solution^^®  and  the  fully 
developed  rotational  flow  solution.  A comparison 
with  the  data  of  Gardon  and  Akfirat,  assuming  that 
the  half  width  is  half  the  jet  width  (i.e. , core  flow 
is  present),  shows  that  the  measured  ground  plane 
pressures  are  between  the  eomputed  developing  flow 
cases  o = . I and  o = .25.  Qualitatively  similar 
results  are  found  in  the  axisymmetric  case  when 
datal2,19  are  compared  with  the  potential  flow  solution 
of  Strand. 20  From  Bower  et  al.'s^  data,  obtained  with 
a curved  confining  plate  around  the  two-dimensional 
jet  e.xit,  the  half  svidth  at  the  impingement  zone  is 
estimated  to  be  1.  IG  jet  widths,  in  reasonable  agree- 
ment with  the  assumption  of  a developing  flow.  Although 
the  ground  plane  pressure  data  differ  from  the  present 
calculations,  a solution  to  the  full  Navicr-Stokes 
equations  by  Bower  et  al."^  shows  no  improvement  (Fig. 
7).  The  confining  plate,  which  can  be  the  source  of  the 
discrepancy  vis  a vis  the  inviscid  calculations,  is  not 
expected  to  influence  the  results  near  the  centerline. 

Fig.  8 demonstrates  that  the  centerline  velocity  decay 
is  closely  represented  by  the  developing  flow  case, 

0 S'  0.5,  while  the  Navier-Stokes  representation"*  still 
produces  poor  results. 

Computation  of  two-dimensional,  oblique  impinge- 
ment produces  the  flow  configurations  depicted  in 
Fig.  9.  As  the  jet  inclination  angle,  0,  decreases 
from  90°  the  stagnation  streamline  shifts,  according 
to  Kqu.ation  (14)  in  order  to  maintain  the  momentum 
Italance.  This  shift  reduces  the  total  pressure  of 
the  stagnation  streamline  (Fig.  10),  although  the 
pressure  distributions  still  appear  to  be  symmetric 
al)out  the  peaks.  The  calculations  show  (ITg.  9)  that 
for  fully  developed  jets,  the  stagnation  streamlines 
approach  the  ground  plane  with  virtually  no  deflection. 

A Taylor  series  expansion  of  the  stream  function  about 
the  stagnation  point  shows  th.it  the  stagnation  line 
incidence  angle,  's  related  to  the  velocity  gradient 

at  the  stiignation  |X)int,  (3u/9x)g,  and  the  jet  vorticity 
distriliution  liy 


tan  0 

s 


2(du/9X) 


(24) 
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Fig.  7 Ground  Plane  Pressure  Distribution, 

Two-Dimensional  Normal  Impingement, 
Fully  Developed  and  Developing  Jet  Profiles 


Fig.  8 Centerline  Velocity  Decay,  Two  Dimensional 
Normal  Impingement,  Fully  Developed, 
Developing  Jet  Profiles 
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Fin.  Coniputalioiis  of  Jot  Half  width  and 
Siannation  Slrctuiilines  (a)  d-  iiU", 
(l»  d 70°,  (c)  (d)  d ao” 


GROUND 
PLANE 
PRESSURE,  p 


-5  -4  -3  -2 


DISTANCE  FROM  ORIGIN,  x 


FiR.  10  Fully  Developed  Flow  Ground  Plane  Pressure 
Distribution  Computations  for  2-D  Oblique 
•let  ImpinRement 


The  data  of  Schauer  and  Eustis^?  Is  compared  with 
the  present  computations  for  the  fully  developed  jet 
by  again  relating  jet  heights  and  half  widths  (Table  4). 
The  observed  stagnation  point  location  (Fig.  11)  is 
in  agreement  with  the  computations  at  0=  70*  but 
deviates  with  reduced  impingement  angles.  Included 
in  the  calculated  results  is  a case  with  3=  20°  in 
order  to  illustrate  that  the  magnitude  of  the  observed 
trends  are  recovered  by  the  computations. 


Table  4 Dimensionless  jet  impingement  zone  height,  z , 
and  free  jet  exit  plane  height,  h , measured  from 
the  ground  plane  along  the  jet  centerline  as  a 
function  of  impingement  angle. 


3 (deg) 

z 

h 

z/h 

90 

.8.  4(i 

12*  2 

,28 

80 

3,51 

12.2 

.28 

70 

3.  57 

12.3 

.29 

(iO 

3,  83 

12.5 

.31 

CiO 

4.  13 

12,8 

. 32 

40 

4,84 

13.5 

• 3(> 

•TO 

5.  78 

14,5 

.40 

_ SCHAUER  8i  EUSTIS'' 

O (STAGNATION  POINT,  x,) 

ROTATIONAL  FLOW 

(PRESENT  RESULTS! 

POTENTIAL  FLOW 


where  ui  (x^)  is  the  vorticity  on  the  stagnation  stream- 
line. This  relationship  which  is  well  satisfied  in  the 
nunierieal  calculations  serves  to  verify  the  computational 
results  and  indicates  that  for  a uniform  (potential)  jet 
the  slagiiation  streamline  is  perpendicular  to  the  ground 
plane.  It  is  instructive  to  compare  further  the  rotational 
and  potential  jet  cases.  The  location,  Xg,  of  the 
resultant  force  is  almost  identical  for  the  two  cases 
over  the  range  '.iO'  <3  < 90°  but  the  deflection  of  the 
uniform  jet  stagnation  streamline  eauses  a signiHcant 
displacement  of  the  maximum  pressure  point22,  .\ 
whereas  Xg  and  Xg  arc  virtually  coincident  for  the  fully 
developed  jet  (Fig.  11).  The  separation,  Xg-Sg,  implies 
a ground  plane  pressure  distributionlO  that  is  asymmetric 
about  the  stagnation  point  (Fig.  12)  in  contrast  to  the 
result  for  the  fully  developed  jet. 


90  80  70  60  50  40  30  20 

JET  INCIDENCE  ANGLE.  0 (DEG) 

Fig.  11  Ground  Plane  Location  of  Itesultant  Force  and 

Stagnation  Point  for  Rotational,  Irrotational 

Jets  (X  , Resultant  Force;  x , Stagnation  Point) 
s s 
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Conclusion 


The  ground  plane  pressure  distribution  data 
exhibit  symmetry  about  the  stagnation  point  as 
suggested  earlier  and  are  well  represented  by  the 
calculations  for  60*  </S < 90' (Fig.  13).  This  is  con- 
sistent with  the  range  over  which  agreement  is  obtained 
for  the  location  of  the  stagnation  point.  Quantitative 
disagreement  between  the  computations  and  data  at 
more  shallow  incidence  angles  is  not  unexpected.  The 
calculated  height  of  the  impingement  zone  is  40%  of  the 
jet  height  at  3 = 30'  (Table  4)  which,  along  with  order 
of  magnitude  arguments  that  can  be  invoked, ®>23 
illustrates  the  increasing  influence  of  the  turbulent 
stress  gradient  terms  and  consequent  breakdown  of 
the  inviscid  rotational  flow  model  at  shallow  impinge- 
ment angles. 


GROUND  PLANE  PRESSURE,  p 


DISTANCE  FROM  STAGNATION  POINT,  x-x, 


Fig.  12  Potential  Flow  Ground  Plane  Pressure 

llistribution  for  2-D  Oblique  Jet  Impingement 


DISTANCE  FROM  STAGNATION  POINT,  (x-x,)/h 


Fig.  13  Comparison  of  Computed  and  Observed  Ground 
Plane  Pressure  Hislributions  for  2-1)  Oblique 
Impingemenl 


Numerical  computations  for  a variety  of 
rotational  jet  impingement  problems  are  given 
here.  From  the  results  and  discussion  it  may  be 
concluded  that: 

(1)  the  inviscid  rotational  flow  formulation  provides 

an  excellent  description  of  the  ground  plane  pressure 
distribution  for  normally  impinging  fully  developed 
a.\isymmetric  and  two-dimensional  jets, 

(2)  ground  plane  pressure  distributions  and  centerline 
velocitj’  decays  for  normally  impinging  developing 
two-dimensional  jets  may  also  be  represented  by 
this  technique  and  show  improvement  with  respect 
to  existing  Navier-Stokes  type  computations, 

(3)  the  structure  of  obliquely  impinging  fully  developed 
two-dimensional  jets  is  such  that  the  stagnation 
streamline  exhibits  little  deflection  and  the  ground 
plane  pressure  distributions  are  symmetric  about 
the  stagnation  point, 

(4)  quantitative  agreement  with  observations  is  obtained 
for  stagnation  point  location  and  ground  plane 
pressure  distribution  for  oblique  jet  impingement 
angles  (3  >G0',  ar.d 

(5)  the  e.\tent  of  the  impingement  zone  increases  as  the 
jet  inclination  angle  drops  (e.g. , at  0 ■-  30'  it  is 
40%  of  the  jet  height)  providing  one  indication  that 
turbulent  stress  gradient  terms  tecomc  increasingly 
relevant  in' the  impingement  zone  for  3<G0'. 
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